Capacities of linear quantum optical systems 
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A wide variety of communication channels employ the quantized electromagnetic field to convey information. 
Their communication capacity crucially depends on losses associated to spatial characteristics of the channel 
such as diffraction and antenna design. Here we focus on the communication via a finite pupil, showing that 
diffraction is formally described as a memory channel. By exploiting this equivalence we then compute the 
communication capacity of an optical refocusing system, modeled as a converging lens. Even though loss of 
information originates from the finite pupil of the lens, we show that the presence of the refocusing system 
can substantially enhance the communication capacity. We mainly concentrate on communication of classical 
information, the extension to quantum information being straightforward. 
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I. INTRODUCTION 

The most prominent candidate for implementing long dis- 
tance quantum communication is undoubtedly represented by 
the electromagnetic field (EMF) |Q]]. Although quantum in- 
formation theory is more commonly represented in terms of 
discrete variables (e.g., qubits), information is most naturally 
encoded in the EMF by means of continuous variables, which, 
in the quantum domain, are described by bosonic degrees of 
freedom. Moreover, all the fundamental quantum informa- 
tion tools and protocols have been demonstrated for continu- 
ous variable systems JHHh, from quantum computationjH] to 
quantum error correction Ipjlal, quantum teleportation 0] and 
quantum key distribution IISl— TlOfl - Here we consider the prob- 
lem of optical quantum communication ifTlULlll . and compute 
the communication capacity, i.e., the maximum rate at which 
information can be reliably transmitted. Although we explic- 
itly consider communication of classical information (3, our 
results are immediately extensible to the case of quantum in- 
formation lfl5ll . 

The most general and simple, although physically rele- 
vant, mathematical model of optical communication line is 
the Gaussian channel |fl6ll , which describes the linear propa- 
gation of the EMF. In the classical domain, the ultimate limits 
for communication via Gaussian channels were provided by 
the seminal work of Claude E. Shannon lfl7tl . In the quan- 
tum domain, the structure of Gaussian channels is notably 
rich IU8I1 . with non-trivial properties in terms of degradabil- 
ity Il9ll and security Ir20l l2lll . However, a full information 
theoretical characterization has been presently achieved only 
for certain families of channels, such as the lossy channel 1221 - 
l24ll . These results have been applied to compute the maxi- 
mum rate of reliable communication via attenuating media, 
as o ptic al fibers, wave-guides, and via free-space propaga- 
tion Il3l I25I I2r3l . Moving along this line we shall provide the 
information theoretical description of the effects of the signal 
propagation through lossy communication channels with lin- 
ear characteristic. After introducing the general methods, we 
consider the example of an optical refocusing system with fi- 



nite pupil, which is schematized as a thin lens which is placed 
between the sender of the message and the receiver under fo- 
cusing conditions. Notwithstanding its relatively simple struc- 
ture this setup captures the basics features of all those situa- 
tions in which a transmitted signal is either focused on a de- 
tector by a suitable optical system prior to the information de- 
coding process, or where it has to be refocused by a suitable 
repeater to allow long distance communication, e. g., b y means 
of parabolic antenna for satellite communication 112711 . 

We explictly discuss how the signal diffraction through the 
optical system can be formally described as a quantum chan- 
nel with correlated noise (memory channel) 12811 . where cor- 
relations are quantified in terms of the associated Rayleigh 
length. Therefore the information theoretical characterization 
of the resulting quantum channel is carried out using tools and 
methods that have been previously applied to quantum mem- 
ory channels 12911 . In this framework the main effect of the sig- 
nal propagation through the optical system is to introduce the 
diffraction of light caused by its finite pupil, leading to band- 
width limited communication lT30tl . However, when compared 
with the free-space communication scheme lT26ll . the presence 
of the refocusing apparatus may yield an improvement in the 
channel capacity of the system. This possibility has been put 
forward in l3ltl . Here we provide a detailed derivation of the 
channel model together with the analysis of results in several 
configurations. Furthermore, we also extend the approach to 
encompass the case of non-monochromatic light. 

The article proceeds as follows. In Sec. [Til we provide the 
general model for communication via lossy multimode quan- 
tum optical channels. In Sec. [TTT1 we introduce a simplest ex- 
ample of such a system, communication through a converg- 
ing lens with finite pupil, and show that it can be formally 
described as a channel with correlated noise. In Sec. [TV] we 
quantize this system, whose communication capacity is com- 
puted in Sec.[V]for the case of monochromatic light. We then 
compare the communication via the refocusing system with 
the free-space setting in Sec.[vT] Section lVTil is devoted to the 
extension to the non-monochromatic case, and Sec. IVIIll is for 
final remarks. 
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H. COMMUNICATION CAPACITY OF QUANTIZED 
LINEAR OPTICAL SYSTEMS 

Consider a linear optical system with a set of M transmit- 
ter modes, labelled by i, and N receiver modes, labelled by j. 
In the case of radio-frequency or microwave communication, 
for example, the transmitter and receiver could be antennae. 
For optical communication, the transmitter could be a laser 
coupled to a telescope, and the receiver could be a telescope 
coupled to a charge-coupled device (CCD) array. Transmitter 
and receiver modes typically have both spatial characteristics 
determined by the optical properties of the setup, and temporal 
characteristics, determined by the frequency and bandwidth of 
the radiation employed in the communication. The transmit- 
tivity matrix Tji gives the fraction of light from the i-th trans- 
mitter mode that is received at the j-th receiver mode. We 
would like to determine the maximum amount of information 
that can be sent for fixed total input power. 

First, consider the purely lossy case, in which noise from 
the environment is neglible. This is the case, for example, for 
free-space optical communication in a thermal background. 
The addition of noise will be considered below. When there is 
just one transmitter mode and one receiver mode, the channel 
is simply the lossy bosonic channel, whose classical capacity 
is known 1I22II : if the loss is r\ and the single-use average pho- 
ton number is h, then the number of bits that can be sent down 
the channel is g(r)n), where 



(x + 1) log 2 (x + 1) — x\og 2 x for x > . 







for x ^ . 



(1) 



The capacity is attained by sending coherent states down the 
channel. If there are M parallel channels, with average photon 
number hi for the i-th channel, then the capacity is simply the 
sum of the capacities for the individual channels. 

In our case, we have a single multimode lossy channel with 
transmittivity matrix Tji, which mixes the input modes to- 
gether. This channel can be transformed into a set of parallel 
channels by using the singular value decomposition. The sin- 
gular value decomposition states that any N x M matrix T 
can be written as 



T = VYU , 



(2) 



where V is an N x N unitary matrix, S is an N x M matrix 
with entries only on the diagonal, and U is an M x M unitary 
matrix. In our case, we can write the transmittivity matrix T 
in components as 



Tji — ^ Vjk^/rjkUki , 



(3) 



where {y/rfk} are the singular values of the transmittivity ma- 
trix. The singular value decomposition shows that any mul- 
timode lossy channel can be decomposed into parallel lossy 
channels with input modes corresponding to the rows of U, 
output modes corresponding to the columns of V, and loss 
factors corresponding to the singular values r]k- 

The singular value input modes can now be quantized using 



annihilation and creation operators ctj , a 



t. 



Xj , Uj 



} = s. 



Similarly, the output modes can be quantized using annihi- 
lation and creation operators bj,bj. [bj,bj,] — Sj,f. To 
preserve the canonical quantization relationships, each input- 
output pair is coupled to a environment mode with annihila- 
tion and creation operators 



(4) 



We see that the singular value decomposition of the multi- 
mode lossy quantum channel renders the channel completely 
equivalent to a set of parallel lossy quantum channels with 
loss factors rjj. The communication capacity of the channel 
when hj photons are transmitted down the i-th singular value 
mode is simply J^j diVj^j)- 

If the j-th singular value mode has average energy per pho- 
ton hQj, then the capacity of the channel with total energy E 
per use is obtained by solving the constrained maximization 
problem with Lagrangian 



yielding the solution 



Vj 



(5) 



(6) 



Here, the Lagrange multiplier fx is chosen to give the 
proper total energy E. In the following, when dealing with 
monochromatic light, we express the energy in terms of num- 
ber of photons, h~Ej (huj). 

Thus, the singular value decomposition allows us to trans- 
form any linear, lossy, multimode channel with transmittivity 
matrix Tji into a set of parallel lossy channels. Quantization 
then yields the capacity of this set of parallel channels. We 
now apply this result to various simple optical settings. In 
particular, we show that the resulting capacity can be signif- 
icantly larger than the capacity which is achieved by using 
'naive' coding and decoding techniques. 



III. THE OPTICAL SYSTEM 

We consider the propagation of light along an optical axis z . 
The input and output signals are identified by the transverse 
light fields at two planes orthogonal to the optical axis, re- 
spectively identified as the object plane and image plane, see 
Fig. QI a )- In classical scalar optics, the monochromatic opti- 
cal fields at the object and image plane are described by scalar 
functions U (r ), Ui(r.i), where r Q and r,; are the cartesian 
coordinates at the object and image plane. For linear systems, 
including the free-space propagation of light, the input/output 
relations at the object and image planes, respectively, are de- 
scribed by a transfer function T(Yi, r D ), such that 



Ui{Ti) 



d 2 Y T{Vi,Y )U (Y ) 



(7) 



3,0 



The quantum version of the relations (0 can be derived by ap- 
plying the canonical quantization procedure and introducing a 
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(b) 



(c) 
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FIG. 1 : (Color online.) (a) Optical communication through an opti- 
cal refocusing apparatus, modeled as a thin lens of radius R and focal 
lenght /. (b) The free-space propagation scenario, (c) An alternative 
scenario in which the lens is replaced by a hole of the same size in 
the absorbing screen. r and r, denote the Cartesian coordinates on 
the object and image planes, respectively. 



proper set of normal modes for the input and output fields (as 
it is done in lfl3ll for the free-space propagation). 

As mentioned in the introduction we model the optical re- 
focusing system as a converging lens of focal length /, lo- 
cated at distance D a from the object plane. Working in the 
thin lens, paraxial approximation, and neglecting aberrations, 
light is collected at the image plane located at distance Di 
from the optical system, where 1 /D a + 1 /Di = 1/ /. Eventu- 
ally the image is magnified by a factor M = £),•/ D Q . Diffrac- 
tion of light is responsible for image blurring and causes loss 
of information. It can be described by introducing an effec- 
tive entrance pupil characterizing the optical system. Denot- 
ing .P(r) the characteristic function of the pupil that encircles 
the lens, the transfer function for the monochromatic field of 
wavelength A is obtained by Fourier transforming l33"ll : 



T(r t 



X 2 D a Di 



d 2 rP(r) e 



(8) 



(throughout the paper l denotes the imaginary unit). In writ- 
ing the above expression we are implicitly assuming that the 
light rays which do not hit the pupil will not reach the im- 
age plane (either because they are scattered out, or because 
they are absorbed by some medium). Such hypothesis is not 
fundamental and could be dropped by adding an extra term 
to Eq. © which accounts for the diffraction of the light rays 
that miss the pupil. Eventually, we notice that, in the paraxial 
approximation, the acquired phase results 



2nD n 



XD n 



,| 2 + | ri | 2 /M) + ^(l + M) . (9) 



For instance, in the case of a circular lens of radius R, 
pupil function is 



the 



P(r) 



1 for 
for 



\t\<R, 
\r\>R, 



and the transfer function reads 

e "?(r s , 



T(r,,r ) 



>i? 2 Jl(27Ti?p) 



X 2 D a D t 



Rp 



(10) 



(11) 



where J\ indicates the Bessel function of first kind and order 
one, and p := |r< — Mr \/(XDi). 



A. Light diffraction as a memory channel 

We shall see that the effects of diffraction on light propaga- 
tion can be described as memory effects in the communication 
channel. Let us recall that a memoryless channel is such that 
its action on different channel inputs are identical and inde- 
pendent. Conversely, the actions of a memory channel, also 
called a channel with correlated noise, at different uses are 
not independent and/or not identical. 

Let us assume that information is encoded at the object 
plane by an array of pixels located at positions r Q (k), with 
the integer k labeling the pixels. Different signals are emitted 
from different pixel positions, which play the role of an array 
of independent channel inputs. It follows from Eq. (0 that the 
input from the fc-th pixel is mapped to an output field at the 
image plane with spatial amplitudes 



Ul k \v i )=T{v h v {k))U^ 



(12) 



Let us notice that, even though the action of the channel is 
identical for all the input pixels, the output fields are not mu- 
tually independent. In fact, a pair of output fields, correspond- 
ing to the the fc-th and fc'-th inputs have a non vanishing spatial 
overlap 



r,[/f>(r,)C/r ; >») 

^(fc) £,(*')* 



(A-') 



that is, 



Ck,k' 



T(T i ,T o (k))T*{T i ,T (k')). 



(13) 



(14) 



The overlap between output signals may cause interference, 
which in turn produces distortion and loss of information in 
the communication via the optical channel. In particular, the 
overlap with the output field generated by the fc-th input sig- 
nal induces noise in the detection of the fc'-th output. If the 
overlap between the two output fields is not negligible — i.e., 
diffraction produces sensible effects — the noise affecting the 
fc'-th output field turns out to be highly correlated with the 
input signal at the fc-th pixel. 

Introducing the dimensionless variable f := r/R (here R 
is the linear extension of the entrance pupil) and using the 
expression in Eq. ^ for the transfer function, we get the fol- 
lowing expression for the output signal overlap 



C 



fc.fe' 



d 2 iP(Rv)e 



(15) 



from which it is evident that the spatial overlap between the 
two output fields is determined by the Rayleigh length of the 
apparatus 



x R 



XDg 
R 



(16) 



4 



For instance, in the case of a circular pupil of radius R, we 
obtain 

Ckk = e M M2n\v (k)~v (k')\/x R ) 
x R \r a (k) - r (k')\ 

In conclusion the Rayleigh length, which quantifies the 
amount of diffraction in the optical system, also quantifies the 
degree of correlations in the optical communication channel. 
The corresponding quantum memory channel manifests inter- 
symbol interference effects, and is qualitatively analogous to 
those studied in JHHl]. 



IV. THE QUANTUM CHANNEL 

The light field at the object and image planes can be quan- 
tized according to standard canonical quantization. In order 
to derive the quantum version of the input/output relations in 
Eq. (0 we first identify a proper set of normal modes at the 
input and output field and proceed along the lines detailed in 
Sec. HI] 

We assume that information is encoded in the object plane 
on a square of length L, creating an image on the image plane 
which (in the geometric optics approximation) is contained in 
a square of size ML. We hence introduce the field variables 

U (n ) := yj<pT e-H^-®*-*)u (T ), 

(18) 

where the integral over r is restricted to the surface of area 
L x L which encircles the object, while the integral over is 
restricted to the surface of area ML x ML which defines the 
receiving screen, and where n Q and are vectors having two 
integer components. The functions U (n ) and C/i(n,) ex- 
press the field components of transverse momentum 2im / L 
and 27rrii/ (ML), respectively at the object and image plane. 
Substituting Eqs. ( fT8l into (0, and using ®, we write the 
input/output relations in the form 

Ui{m) = T *H,no U (n ) , (19) 
where the transfer matrix T nino reads 

^ = i// A ^ r) ' (20) 

with 

x JJd^e- 12 ^^^) ^ . (21) 



Finally, we consider the singular value decomposition of the 
transfer matrix, 

T ni ,n = ^ni,n ^/Vn^n^a ! (22) 

n 

where, as in Eq. (|2), U nyIlo , V ni . n are unitary matrices, and 
{V^/n} are tne singular values of the matrix T nii „ o , taking 
values in the interval [0, 1]. A set of input and output field 
variables are hence defined as follows: 

U a (n) := Y, M n,o„ U(n ) , (23) 
C7,(n) := ^V^E^m). (24) 

It follows from Eq. ( l22b that they satisfy the identities 

Ui(n) = V^o(n) • (25) 

In other words, the field variables U (n) are independently, 
but non necessarily identically, transmitted to the output vari- 
ables f/j(n). The effect of the channel is to attenuate the n-th 
variable by a factor r/ n . 

We can now promote the output and input field variables to 
the rank of quantum operators, by substituting 

U (n) -> an, U:{n)^ a+ , (26) 

Ui(n) -> ^/KJj2b a , U*(n) -> ^/KJj2bl , (27) 

where w = 2irc/A is the frequency, and imposing the canoni- 
cal commutation relations: 

[a n , a£,] = S^ n , , (28) 
K,bl] =«„,„/. (29) 

The preservation of canonical commutation relations requires 
to invoke a set of canonical noise variables {£n,£ n } and to 
write the quantum version of Eq. ( |25T > as 

K = 07n On + \J 1 - ?7n ■ (30) 

This set of input/output relations, together with their hermitian 
conjugates, characterizes, upon evaluation of the transmissiv- 
ities rj n , the quantum description of the optical channel. 



V. CAPACITY OF THE OPTICAL QUANTUM 
COMMUNICATION CHANNEL 

The input/output relations for the quantum description of 
the channel which have been derived in Eq. d30b formally 
qualify the optical system as a broadband lossy channel: a 
multimode channel in which a collection of bosonic modes is 
transmitted with corresponding efficiencies. This model has 
been characterized from the information-theoretical viewpoint 
in l22il23Tl . where the capacities of the channel for transmit- 
ting classical and quantum information have been established. 
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The number of modes which are transmitted through the 
channel is virtually infinite. Actually, due to the finiteness of 
the entrance pupil, only a finite number of modes has nonzero 
transmissivity. The values of the transmissivities r/ n can be 
numerically estimated by first computing the elements of the 
transfer matrix T ni . no , and then its singular values. Indeed, 
the matrix in Eq. (fJTJ can be rewritten in terms of the ratio 
L/xr, and an analytical solution for the values of the effec- 
tive transmissivities can be deduced in the following limits, 
namely 



L < x R 



L > Xr, 



FAR FIELD, 



NEAR FIELD. 



(31) 



(32) 



First of all, we notice that in the far-field case Eq. (I2TI 1 
reads: 



A n ,,n ~ ML 2 <5 nij0 <5n o ,0 : 



(33) 



where S denotes the Kronecker symbol. This expression in 
turn leads to 



(34) 



that is, in the far-field regime only one mode is transmitted 
through the optical channel, and it is attenuated by a factor 
~ 7r 2 (L I xr) 4 . At least for optical frequencies it is reasonable 
to assume that the noise modes in Eq. (l30T > are not populated. 
Under this assumption, using the result of 112211 . we are led 
to the following expression for the classical capacity of the 
optical channel: 



C s = g 



n 2 L 4 



(35) 



where h denotes the mean number photons at the object plane. 
The above calculations refer to the case of scalar field. Po- 
larization can be included by doubling the total number of 
modes. In this case Eq. ( f35l > will then be replaced by 



2<? 



n 2 L 4 
2x R 



(36) 



where the factor 2 which multiplies the g-function comes from 
the fact that now there are two modes which can efficiently 
propagate, while the extra factor 1/2 in the argument of g 
comes from the fact that the available energy must be equipar- 
titioned between them. Analogously, the quantum capacity 
can be computed according to 12311 . 

Let us now consider the near-field limit, L ^ xr. In this 
regime we can approximate 



An„„ (r) ~ X 2 D Di6 ni , n J^ (r 



L 



(37) 



where 5^ indicates the two-dimensional Dirac function, 
follows that 



(38) 



Within this approximation the transfer matrix is diagonal, 
hence the singular values coincide with its diagonal entries. 
For a circular pupil of radius R, we obtain the following val- 
ues for the transmittivities: 



1 for |n| < L/xr . 
for jnj > L/xr 



(39) 



We deduce that, in the near-field regime, the number of trans- 
mitted modes per surface unit is approximatively equal to 
7r/x|, each being transmitted with approximatively unit ef- 
ficiency. Equivalently we can say that the total number of 
modes which are transferred with unit efficiency is equal to 



v ~ ir(L/xRf 



(40) 



We are now in the condition of computing the capacity. De- 
noting by n the average number of photons impinging on the 
surface, the capacity has the form 13411 



C n f = v g(n/ v) = g n 



(41) 



where we used the fact that the maximum transfer is obtained 
when n is equipartitioned among all transmitted modes. It is 
worth stressing that, since the modes employed in the trans- 
mission are perfectly transmitted, the classical capacity (ex- 
pressed in bits) coicides with the quantum capacity (expressed 
in qubits). As before, our result can be generalized to include 
also the polarization degree of freedom obtaining 



2-kL? ( x 

g ' 



x R " \2ttL 2 



(42) 



To evaluate the capacity for a generic value of the ratio 
L/xr, it is necessary to numerically diagonalize the transfer 
matrix T. For the sake of simplicity, here we present an ex- 
ample of lower dimensionality, in which information is en- 
coded at the object plane along an infinite strip of size L, and 
diffraction is caused by an infinitely long slit of size 2R. The 
analysis of such ID system is analogous to the 2D one, yield- 
ing the following expressions for the far-field and near-field 
capacities. In the far-field limit we have 



(ID) 



and in the near-field limit 



'2L 

g\ — n 

XR 



c. 



(ID) 
nf 



2L (Xr 

— g ' — 

Xr 



\2L 



(43) 



(44) 



Figure [2] shows the transmissivities for several values of the 
ratio L/xr for the one-dimensional problem. Figure [3] shows 
the capacity as function of the ratio L/xr, compared with the 
limiting expressions of Eqs. ( |43l and (l44l i. 



VI. ENHANCED QUANTUM COMMUNICATION VIA 
OPTICAL REFOCUSING 



For a fair comparison between the optical communication 
through the optical system and the free-space one, we consider 
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FIG. 2: The effective transmissivities rj n , for different value of the 
ratio L/xr. Dots: L/xr = 0.1 (far-field). Crosses: L/xr — 1. Cir- 
cles: L/xr — 10 (near-field). The figure refers to a one-dimensional 
setting, in which information is encoded at the object plane along an 
infinite strip of size L, and diffraction is caused by an infinitely long 
slit of size 27?. 




FIG. 3: The plot shows the capacity (in the monochromatic case) as 
function of the ratio L/xr, for n — 4. The figure refers to a one- 
dimensional setting, in which information is encoded at the object 
plane along a straight line of length L, and diffraction is caused by 
an infinitely long slit of size 2R. The solid line is the exact value 
of the capacity computed by numerical evaluation of the set of effec- 
tive transmissivities. The dashed line indicates the approximation for 
near-field limit (L/xr ^> 1). The dotted line is the approximation 
for the far-field limit (L/xr <C 1). The inset is a magnification of 
the far-field region. 



the case of free-space communication under the hypotheses 
that light is emitted by an object of surface L 2 , propagates by 
a distance D = D a + Di = D (l + Ad), and is finally de- 
tected on a surface of size (ML) 2 . This setting is depicted 
in Fig. |TJb). The light propagation, both in the classical and 
quantum regimes, is characterized by the Fresnel number as- 
sociated to this setting, i.e., T = AI 2 L i /(XD) 2 (see JH] 
and references therein). In the far-field limit, J « 1, only 
one mode is transmitted, with a corresponding transmissiv- 
ity equal to the Fresnel number. On the other hand, in the 
near-field limit, T ^ 1, T equals the number of modes 
which are transmitted with unit transmissivity. For these two 
regimes we can hence compute the free-space classical capac- 
ity according to l22tl . By comparison with Eqs. ([35), fiTK 



it follows that the presence of the radius-i? lens enhances 
the classical capacity over the free-space propagation only 
for not too large M, i.e., for not too large imaging plane 
(specifically when ttR/xr > M/(l + AI) in the far-field, 
and when tt^ 2 R/L > M/(l + M) in the near- field limit). 
This rather counterintuitive effect originates from the simpli- 
fying assumption we made in writing Eq. ([8]) which implicitly 
states that the only light rays reaching the imagine plane are 
those which passes through the pupil, while the other are lost. 
Clearly the presence of such loss mechanism (which is not ac- 
counted for in the free-space calculation) is uninfiuential as 
long as the image plane is small, while it becomes relevant 
for large imaging screen: this is exactly where the free-space 
starts to outperform the propagation through the lens. 

Below we present a comparison between the performances 
of the refocusing apparatus, scenario (a), and those of the free- 
space propagation case, scenario (b). The goal is to produce 
results which are not affected by the approximation we made 
in writing Eq. ([8]l (i.e., the fact that we have implicitly as- 
sumed that all photon which do not hint the lens will not be 
transferred on the imaging plane). In order to do that, we also 
consider a third scenario, denoted by (c), which is depicted in 
Fig. He). 



A. Far-field regime for scenario (a) 

Let us first consider the case in which the scenario (a) is op- 
erated in the far-field regime d3l~l ), which according to Eq. ( fT6l ) 
corresponds to have 



A» 



LR 
~D~o 



(45) 



As already seen in the previous section, under this condition 
the scenario (a) is characterized by having a single mode trans- 
mitted from the object plane to the image plane with an effec- 
tive transmissivity 



V(a) = tt 2 (L/x r ) 4 



( LR 
\XD~o 



< 1 . 



(46) 



In the scenario (b), the field propagates freely from the object 
plane to the image plane. In this case the far-field regime is 
equivalent to impose 



LML L 2 



M 



D 



D n M+l 



(47) 



which in principle is independent from the far-field condi- 
tion d45l ) for (a). Under condition d47| i also (b) will admit 
the propagation of a single mode which is now attenuated by 



L Li 

Id 



AI 



AI + 1 



v- 



XD 



< 1 , (48) 



where L a = L is the dimension of the object plane, Li = 
AIL the dimension of the image plane, and D = D a + Di 
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is the distance between them. Then, assuming that both (l45l l 
and d47T > hold, the ratio between the transmissivities is 



( R2 



\XD 



M + l 
M 



(49) 



As already noticed this is not always larger than one: indeed 
there is a condition which M, R, D a , A have to satisfy for this 
to happen. We now show that such condition is equivalent 
to impose that the loss induced by the pupil should be negligi- 
ble. To do so consider the alternative scenario (c), in which the 
field propagate from the object plane to the image plane while 
instead of the lens we have just a hole in the absorbing screen 
of the pupil. This configuration can be treated by analyzing 
separately the free-space propagation from the object plane 
to the absorbing screen (o — > s) and the free-space propaga- 
tion from the absorbing screen to the imaging plane (s — > i). 
We notice that the condition Eq. d45b imposes that both these 
propagations take place in the far-field regime lf35ll . Therefore 
we can conclude that, under the condition ((45), in the scenario 
(c) there will be a single propagating mode (o -+ s -+ i). It is 
attenuated by a factor 



M _ 



in the section (o — > s) and by 



LR 



"(c) 



( RML\ 



LR 
ATX 



(50) 



(51) 



in the section (s — > i). The overall attenuation cannot be 
bigger than the product of the two factors, i.e., 



(o- 



•s) (s- 



LR 
\D~o 



(52) 



Now, if the presence of the absorbing screen around the pupil 
has to be negligible, then we must have that the loss that it 
induces are equal to the one of free-space propagation. Since 
Vic) ^ V(a), this implies that the regime in which we can ne- 
glect the effect of absorbing screen is exactly the one in which 
n is greater than 1. In other words the detrimental effects 
we see for r\ < 1 simply correspond to the screen's absorp- 
tion. As a final remark we also observe that the condition 
r\ > 1 plus the far-field regime (05J for (a) enforce the far- 
field regime d47| i for the free-space propagation, indeed, 



XD M + l 



L 2 M 




> 1 (53) 



We now compare the performances of scenario (a) and (b) 
in terms of capacities by computing the ratio 



Gi 



C(a) _ g(V(g)n) _ g(rir/ {b) n) 



(54) 



This is a monotonic function of n, with G\ ~ r\ for n <C 1 
(faint signal limit) and G\ ~ 1 for n 3> 1 (semiclassical 
limit). The enhancement in the transmission rate provided by 
the optical refocusing system persists in the presence of back- 
ground thermal noise. In such a case, by encoding classical 
information into coherent states, the expression g{rjfi) has to 
be replaced by g(rjn + nth) ~ <?(^th), where nth denotes the 
mean number of thermal background photons per transmitted 
mode Ql, yielding d ~ r\ even in the very noisy limit 
n th > fi. 



B. Near-Field regime for scenario (a) 

The near-field regime for the scenario (a) is defined by the 
inequality ( f32t which rewrites also as 



LR 



(55) 



From the discussion of Sec. [V] we know that in this regime 
the scenario (a) is characterized by having a collection of 
modes which are perfectly transmitted from the object plane 
to the image plane, 



'(a) 



= k(L/x R ) 2 = 7T 



LR 
ATT 



> 1 



(56) 



In a similar way the near-field regime for the scenario (b) takes 
place when 



A< 



L n Li 



D 



L 2 M 
IX M + l ' 



(57) 



is 



independent from the near-field condition for 

modes 



which 

(a) (155b . Under this condition also (b) admits vn,\ 
which are perfectly transferred, with 



V{b) 



L Li 



M 



M + l 



XD 



> 1 . (58) 



Assume then that both ( l55l l and ( 1571 ) hold and define the qual- 
ity factor 



''2 



(M + l 



M = 

v {b) \ M 



(59) 



As in the case of Sec. lVI Al if there are no losses introduced 
by the absorption of the rays propagating outside the lens, we 
must have T2 ^ 1. To see this let us consider what happens in 
the scenario (c). We first notice that the condition (|55l ) guar- 
antees that both the propagations (o -+ s) and (s — > i) are in 
the near-field regime. The number of modes that they allow 
for perfect propagation is given by 



(o— >s) (s 
(c) (c 



LR 
XL\, 



>) 



(60) 



Notice that they are identical, due to the fact that L; / 1 Di = 
L/D Q . This implies that V( c ) ^ v^ a y Then, it is clear that 
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the presence of the pupil is negligible only when vr c \ is larger 
than vn,), i.e., r 2 ^ 1, therefor proving the thesis. 

We remark that the condition r-z ^ 1, together with the 
near-field condition for the scenario (a) [Eq. ( f55l ll is not suffi- 
cient to garantee the near-field condition for the scenario (b). 
Indeed we notice that 



M 



LR 



XD M + 1 \D KM + 1R 



M L 



(61) 



which by itself does not imply the near-field condition for (b), 
i.e., Eq. ( 1571 ). That yields the possibility to have (a) in near- 
field and (b) in either near-field or far-field. 
Let us examine the two cases. 

• First assume that both (a) and (b) are in the near-field 
regime (of course under the constraint that r 2 > 1 to 
exclude the absorption by the pupil). The ratio between 
the capacities becomes 



G, 



C(a) _ V(g) g(n/v( a )) 

C(b) V(b) g{n/v( b )) 



r-2 



ffOW^a)) ' 



(62) 



which is a monotonic function of n, with G 2 — 1 for 
n <C 1 (faint signal limit) and G 2 — ?"2 for n ^> 1 
(semiclassical limit). 

• Then assume that (a) is near-field and (b) is far-field, 
i.e., 



L 2 M LR 

^ji7TT <<A<< ^ 



(63) 



Under these conditions one can immediately verify that 
the presence of the absorbing screen does not affect the 
performances of the scenario (b), and we can make a 
fair comparison. In this case the gain becomes 



G, = 



G(a) _ V{a) g{n/V(g)) 



(6) 



007(6)' 



(64) 



VII. 



giving G3 ~ l/?7(b) 3> 1 for n <C 1 (faint signal limit) 
and G3 ~ Vt a \ S> 1 for b > 1 (semiclassical limit). 
In the presence of noisy thermal environment lfl6ll . the 
expressions g(i]fi) and g(njv) have to be replaced by 
g{r\n + n t h) - g{fhh) and g(njv + n th ) - g{n t h), re- 
spectively. The advantages of optical refocusing hence 
persist even in the very noisy limit nth ^ n, yielding 
G 3 ~ l/??(6). 



COMMUNICATION WITH NON MONOCHROMATIC 
LIGHT 



Till now we have considered the case of monochromatic 
light, we now want to compute the communication capacity of 
the optical system by assuming non-monochromatic light. An 



input signal over a time interval T can be expanded in terms 
of monochromatic components at frequencies u>j = 2irj/T, 
with j = 0, 1, ... 00. If the optical system is characterized by 
a certain bandwidth extending from to £1 + 50,, only a finite 
number of components are transmitted, corresponding to the 
frequencies us n such that il ^ ljj ^ il + 60.. Each monochro- 
matic components contributes with a term as in Eq. (1411) . with 
a-R = XjD /R = 2ircD / (ojjR). In the following we will 
work under the assumption that the frequency modes are ei- 
ther all in the far-field regime or all in the near-field one. 



A. Far-field 

If all the transmitted frequency modes fulfill the far-field 
condition Pil l, we must have 

L < x R = 2ircD /(uj j R) , (65) 

for all LOj . Consequently, the results of Eq. ( 1351 ) can be applied 
to the whole spectrum. This allows us to derive the following 
expression for the capacity 

Gff = 2J 9 i (66) 



where 



LR 



\2itcD 



(67) 



and the parameter fij counts the average number of photons at 
frequency uij . If a mean power P is employed, the parameters 
hj ought to obey the constraint 



rp 



huij rij 



P. 



(68) 



The maximization over photon number distributions satisfy- 
ing the input energy constraint can be worked out by Lagrange 
method, yielding the optimal photon number distribution 



aw; 



7^ 



1 



(69) 



with fi being the Lagrange multiplier. For sufficiently large 
T we can approximate the summations with integrals, and the 
channel capacity reads 



T 



n+sn 



dw g 



2^6* - 1 



(70) 



where the value of the Lagrange multiplier is determined by 
the implicit equation 

r n+sn ^ 

,r5 



P 



h 

2(XK 



1 



(71) 



A closed form for the classical capacity can be found in the 
narrow-band limit, 50. <C O, in which we obtain the approxi- 
mate expression 



a 



T50 (2TrPa0 3 
~2^~ 9 \ MO 



(72) 
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B. Near-field 

Let us now assume that all the frequency modes fulfill the 
near-field condition (|32| |. i.e., 



L > x R = 2wcD /(uijR) 



(73) 



We can hence apply Eq. (f4TT > to the whole spectrum, which 
allows us to write the channel classical capacity as follows: 



Cnf 



E 



v i9\ — 



(74) 
(75) 



where Uj counts the number of transmitted modes of the fre- 
quency cjj [Eq. d40bl. hj counts the average number of photon 
of that frequency, and where we have introduced 



/3:=tt 



/ LR 
\2ttcD 



(76) 



The optimization over the photon-number distribution under 
the constraint of total power P yields the optimal photon- 
number distribution 



(77) 



with fj, being the Lagrange multiplier. For sufficiently large 
T we can approximate the summations with integrals, and the 
channel capacity reads 



~ 2irq 3 



n+sn 



doj uj 2 g 



1 



dx x 2 g 



qD. 



1 



e x - 1 



, (78) 



where we rescaled the Lagrange multiplier introducing the 
quantity q = In (2)/ift/T, that is determined by the implicit 
equation 



P 



ph /■«("+*") , x 3 



2irq 4 



dx ■ 



qQ 



1 



(79) 



We can single out two limiting situations, the narrowband and 
the broadband limits, for which analytical expressions can be 
obtained. 



1. Narrowband limit 



In the narrowband limit, 8fl <C fi, we obtain 



CI? ~ — T STSSl g 
ni 2n y 



( 2irP 



(80) 



This can be casted in a more familiar form by noticing that in 
this limit the power P can be expressed as 



p4(^) mB(n) 



(81) 



where 8QT/(2n) is the total number of frequencies and n(£Y) 
is the density of mean photon number at frequency 0. Replac- 
ing this into Eq. (f80b and using /30 2 = 7t(L/xr) 2 , where xr 
is the Rayleigh length of the frequency O, we get 



b ^ (50T\ L 2 



(82) 



This expression shows that C£f coincides with the single fre- 
quency capacity d4Tb , multiplied by the total number of fre- 
quencies 5QT/{2it). 



Broadband limit 



In the broadband limit, we set O 



p 



bb 



°nf 



2irq 

/3T 
2ixq c 



sn -> 

4 , 



-oo, and we have 
(83) 
(84) 



with J-(z) and 9(z) being the following decreasing functions, 



Hz) 



dx ■ 



x 



e" — 1 

2 



da; x g 



1 



(85) 
(86) 



An analytical solution can be obtained by approximating qil 
to zero. Notice that this is not formally correct if we want to 
preserve the condition ( f73l for all frequencies of the spectrum: 
however since J-(z) and Q(z) are smooth in the proximity of 
3 = the error becomes negligible small. By close inspection 
of the Eqs. ( f85l > and d86b however, it follows that the approx- 
imations F(z) ~ -7^(0), Q{z) ~ G(0) are justified if z < 1, 
i.e., for qil < 1, which in turn implies: 



(87) 



Under this assumptions it follows that the condition (l73l ) is 
satisfied in the semiclassical regime P/(Ml 2 ) 1. Conse- 
quently we can write 



P 



Cbb 
nf 



Pn 3 h 
30 9 4 
/3tt 3 T 
45g 3 



log 2 e 



which yields 



n 



bb 
nf 



3\/l5 



2ttP 



3/4 



log 2 e 



(88) 
(89) 

(90) 
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These expressions are obtained by noticing that the integrals 
in ( [85t , ( l86l l can be written in terms of the Bernoulli numbers 
Bk by means of the identity (see, e.g., i3) 

f^_ = ( _ ir -.(M 2 "^. (91) 

Jo eP x -l \pj in 

By comparing d90t with the capacity of a multimode Gaussian 
bosonic channel we notice that the scaling in P is 

now changed (there it scales with as P 1 / 2 ). This is due to the 
fact that in our case each frequency has multiple degeneracies. 

VIII. CONCLUSION 

We have computed the capacity of quantum optical commu- 
nication through an optical system characterized by a finite 
entrance pupil. Our calculations provide general bounds on 
the efficiency of quantum optical communication taking into 
account the effects of light diffraction. This models a rather 
general situation in long distance communication, where re- 
peaters play the role of the optical system used to refocus the 
signal. More generally, any transfer of information which em- 
ploy quantum degrees of freedom of light — from quantum 
key distribution ja- Holl to quantum imaging i38ll as well as 
quantum discrimination problems IB9M42I1 and quantum read- 
ing lf43l - l45ll — requires the propagation through an optical 
system, and it is hence limited by diffraction. 

We have argued that, when the optical system is used for 
information transmission, the effects of diffraction can be for- 



mally described as a quantum channel with correlated noise. It 
follows that correlated noise may affect all the quantum infor- 
mation protocols requiring propagation and detection of quan- 
tum light. 

Finally, this formal equivalence has allowed us to apply 
tools that were developed in the framework of quantum mem- 
ory channel characterization to the diffraction problem. In 
particular, this has allowed us to show that, under certain con- 
ditions, a converging optical apparatus can be used to achieve 
higher transmission rates than the free-space field propaga- 
tion. The tradeoff between loss and diffraction determines the 
conditions under which the intuitive benefits of optical refo- 
cusing can be rigorously proven. 
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